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Abstract — Feedback  control  of  piecewise  smooth  discrete¬ 
time  systems  that  undergo  border  collision  bifurcations  is 
considered.  These  bifurcations  occur  when  a  fixed  point  or  a 
periodic  orbit  of  a  piecewise  smooth  system  crosses  or  collides 
with  the  border  between  two  regions  of  smooth  operation 
as  a  system  parameter  is  quasistatically  varied.  The  goal  of 
the  control  effort  in  this  work  is  to  modify  the  bifurcation 
so  that  the  bifurcated  steady  state  is  locally  attracting  and 
locally  unique.  To  achieve  this,  Lyapunov-based  techniques  are 
used.  A  sufficient  condition  for  nonbifurcation  with  persistent 
stability  in  piecewise  smooth  maps  of  dimension  n  that  depend 
on  a  parameter  is  derived.  The  derived  condition  is  in  terms 
of  linear  matrix  inequalities.  This  condition  is  then  used  as  a 
basis  for  the  design  of  feedback  controls  to  eliminate  border 
collision  bifurcations  in  piecewise  smooth  maps  and  to  produce 
desirable  behavior. 

I.  Introduction 

Stabilizing  feedback  control  laws  for  piecewise  smooth 
discrete-time  systems  exhibiting  border  collision  bifurca¬ 
tions  are  developed.  The  class  of  piecewise  smooth  systems 
considered  constitutes  systems  that  are  smooth  everywhere 
except  along  borders  separating  regions  of  smooth  behavior 
where  the  system  is  only  continuous.  Border  collision 
bifurcations  are  bifurcations  that  occur  when  a  fixed  point 
(or  a  periodic  orbit)  of  a  piecewise  smooth  system  crosses 
or  collides  with  the  border  between  two  regions  of  smooth 
operation.  The  term  border  collision  bifurcation  was  coined 
by  Nusse  and  Yorke  [1],  Border  collision  bifurcation  had 
been  studied  in  the  Russian  literature  under  the  name  C- 
bifurcations  by  Feigin  [2],  [3],  The  results  of  Feigin  were 
introduced  to  the  Western  literature  in  [4]. 

Border  collision  bifurcations  (BCBs)  include  bifurcations 
that  are  reminiscent  of  the  classical  bifurcations  in  smooth 
systems  such  as  the  fold  and  period  doubling  bifurcations. 
Despite  such  resemblances,  the  classification  of  border 
collision  bifurcations  is  far  from  complete,  and  certainly 
very  preliminary  in  comparison  to  the  results  available  in  the 
smooth  case.  In  smooth  maps,  a  bifurcation  occurs  from  a 
one-parameter  family  of  fixed  points  when  a  real  eigenvalue 
or  a  complex  conjugate  pair  of  eigenvalues  crosses  the  unit 
circle.  In  piecewise  smooth  (PWS)  maps,  on  the  other  hand, 
a  discontinuous  change  in  the  eigenvalues  of  the  Jacobian 
matrix  evaluated  at  the  fixed  point  (or  at  a  periodic  point) 


occurs  when  the  fixed  point  hits  the  border.  As  a  result,  bor¬ 
der  collision  bifurcations  for  piecewise  smooth  systems  in 
which  the  one-sided  derivatives  on  the  border  are  finite  are 
classified  based  on  the  linearizations  of  the  system  on  both 
sides  of  the  border  at  criticality.  The  phenomenon  of  border 
collision  bifurcations  has  been  observed  both  numerically 
and  experimentally  in  many  systems,  in  applications  such 
as  power  electronic  devices  [5],  [6],  [7]  and  in  studies  of 
grazing  impact  in  mechanical  oscillators  [8],  [9]. 

There  is  little  past  work  on  control  of  BCBs  [10],  [11], 
[12],  The  control  method  of  [10],  [11]  is  based  on  the 
classification  scheme  of  BCBs  that  was  given  by  Feigin  [4], 
Flowever,  since  Feigin  didn’t  give  conditions  for  specific 
scenarios,  the  results  of  [10],  [11]  do  not  address  stabi¬ 
lization.  Moreover,  references  [10],  [11]  use  a  trial  and 
error  approach  that  doesn’t  provide  analytical  conditions  for 
existence  of  controllers.  In  our  recent  work  [12],  feedback 
control  of  BCBs  in  one-  and  two-dimensional  PWS  maps 
has  been  considered. 

In  the  present  paper,  a  sufficient  condition  for  nonbifur¬ 
cation  with  persistent  stability  in  piecewise  smooth  maps 
of  dimension  n  that  depend  on  a  parameter  is  derived. 
That  is,  a  condition  under  which  the  PWS  map  possesses 
a  locally  asymptotically  stable  fixed  point  which  is  also 
the  locally  unique  attractor  for  all  values  of  the  bifurcation 
parameter  in  a  neighborhood  of  the  critical  value.  This 
condition  is  derived  using  Lyapunov-based  methods  and  is 
given  in  terms  of  linear  matrix  inequalities  (LIMs).  The 
derived  condition  is  then  used  as  a  basis  for  the  design  of 
feedback  controls  to  eliminate  BCBs  in  piecewise  smooth 
maps  and  to  produce  desirable  behavior.  The  analysis  and 
control  methodology  presented  in  this  paper  have  recently 
been  applied  to  a  model  of  cardiac  arrhythmia  [13]. 

The  paper  proceeds  as  follows.  In  Sec.  II,  brief  back¬ 
ground  material  on  BCBs  is  given.  In  Sec.  Ill,  Lyapunov- 
based  analysis  of  PWS  maps  undergoing  BCBs  is  presented 
and  a  sufficient  condition  for  nonbifurcation  with  persistent 
stability  is  derived.  In  Sec.  IV,  the  results  of  Sec.  Ill  are 
used  in  the  synthesis  of  stabilizing  feedback  control  laws 
and  numerical  examples  that  demonstrate  the  results  are 
given.  Concluding  remarks  are  collected  in  Sec.  V. 


II.  Background  on  Border  Collision 
Bifurcations 


Consider  the  one -parameter  family  of  piecewise  smooth 
maps 


f(x,p) 


f ),  x&Ra 

fB(x ,£/),  xGRb 


where  / :  R"+1  — >  K"  is  piecewise  smooth  in  x;  /  is  smooth 
in  x  everywhere  except  on  the  border  (hypersurface  T) 
separating  Ra  and  Rb  where  it  is  only  continuous,  /  is 
smooth  in  p  and  R  \ ,  R/>  are  the  two  (nonintersecting)  regions 
of  smooth  behavior.  Of  great  interest  is  the  study  of  the 
dynamics  of  /  at  a  fixed  point  (or  a  periodic  orbit)  near  or 
at  the  border  T.  If  the  fixed  point  (or  periodic  orbit)  is  in  R,\ 
(respectively  Rb)  and  is  away  from  the  border,  then  the  local 
dynamics  is  determined  by  the  single  map  /a  (respectively 
fs).  If,  on  the  other  hand,  the  fixed  point  is  close  to  the 
border,  then  jumps  across  the  border  can  occur  except  in  an 
extremely  small  neighborhood  of  the  fixed  point.  Therefore, 
for  operation  close  to  the  border,  both  /a  and  fs  are  needed 
in  the  study  of  the  possible  behavior.  For  a  fixed  point 
at  or  near  the  border,  the  dynamics  is  determined  by  the 
linearizations  of  the  map  on  both  sides  of  the  border. 

Various  types  of  BCBs  occur  in  (1)  as  the  bifurcation 
parameter  p  is  varied  through  a  critical  value  [1],  [14], 
[15],  [4].  Such  bifurcations  occurring  in  the  map  (1)  can 
be  studied  using  the  piecewise-linearized  representation  [4] 


x{k T  1)  :=  F{x{k),p) 


Ax(k)  +  bp,  x i  ( k )  <  0 
Bx{k)  +  bp,  xi  (k)  >  0 


where  A  is  the  linearization  of  the  PWS  map  /  in  Ra  at  a 
fixed  point  on  the  border  approached  from  points  in  R  \  near 
the  border  and  B  is  the  linearization  of  /  at  a  fixed  point 
on  the  border  approached  from  points  in  Rb  and  b  is  the 
derivative  of  the  map  /  with  respect  to  p.  The  coordinate 
system  is  chosen  such  that  the  sign  of  the  first  component 
of  the  vector  x  determines  whether  x  is  in  Ra  or  Rb  (a 
transformation  to  the  form  (2)  is  given  in  [16]).  If  xi  =  0, 
then  x  is  on  the  border  separating  Ra  and  Rb-  The  continuity 
of  F  at  the  border  implies  that  A  and  B  differ  only  in  their 
first  columns. 

The  classification  of  BCBs  depends  on  the  eigenvalues 
of  A  and  B  [4].  A  complete  classification  of  BCBs  is  only 
available  for  one-dimensional  PWS  maps  [17],  [18].  For 
two  dimensional  PWS  maps,  some  results  are  available  that 
only  address  a  class  of  two-dimensional  PWS  maps  [1], 
[19],  [15],  [20], 

Although  Feigin  [4]  studied  general  n  -dimensional  PWS 
maps  exhibiting  border  collisions,  only  very  general  condi¬ 
tions  for  existence  of  a  fixed  point  and  period-2  solutions 
before  and  after  the  border  were  given.  The  classification 
scheme  of  [4]  does  not  give  any  information  about  stability 
or  uniqueness  of  fixed  points  or  period-2  orbits  involved  in 
the  border  collision  bifurcation  nor  does  it  give  information 
about  higher  period  periodic  orbits  or  chaos  that  might  be 
involved  in  a  border  collision  bifurcation.  Therefore,  in  this 


paper,  one  of  the  main  goals  is  to  develop  a  sufficient 
condition  for  nonbifurcation  with  persistent  stability  that 
can  be  used  in  the  design  of  stabilizing  feedback  control 
laws.  This  is  done  in  the  next  section  using  Lyapunov-based 
techniques. 

III.  Lyapunov-Based  Analysis  of  Piecewise 
Smooth  Maps 


Recently,  many  researchers  have  studied  stability  of 
a  fixed  point  of  switched  discrete-time  linear  systems 
(e.g.,  [21],  [22],  [23])  as  well  as  continuous  time  switched 
systems  (e.g.,  [24]).  In  all  the  referenced  studies,  Lyapunov 
techniques  were  used  to  obtain  sufficient  conditions  for  sta¬ 
bility  of  the  fixed  point  (or  equilibrium  point)  of  a  piecewise 
linear  system.  For  instance,  in  [24],  [22],  quadratic  as  well 
as  piecewise  quadratic  Lyapunov  functions  were  used  in 
the  analysis  of  stability  of  switched  systems  and  also  in  the 
synthesis  of  stabilizing  controls.  The  author  is  unaware  of 
any  previous  study  using  Lyapunov  methods  to  analyze  the 
dynamics  of  switched  systems  depending  on  a  parameter. 
Here,  a  quadratic  Lyapunov  function  is  used  to  study  border 
collision  bifurcations  in  PWS  maps  and  to  obtain  sufficient 
conditions  for  nonbifurcation  with  persistent  stability  in 
such  maps. 

Consider  the  piecewise-linearized  representation  of  PWS 
maps  given  in  (2).  The  sign  of  the  first  component  of  the 
vector  x  determines  whether  x  is  in  R,\  or  in  Rb-  If  xi  =  0, 
then  x  is  on  the  border  separating  Ra  and  R/j.  The  continuity 
of  F  at  the  border  implies  that  A  and  B  differ  only  in  their 
first  columns.  That  is,  ap  =  bjj ,  for  j  ^  1,  where  A  =  [ ciij ] 
and  B  =  [bij\. 

Assume  that  1  ^  a(A),  1  ^  o(B)  (i.e.,  both  I— A,  I  —  B  are 
nonsingular).  Formally  solving  for  the  fixed  points  of  (2) 
yields  xa(/j)  =  (7  —  A)~lbp  and  xs(p)  =  (I  —  B)~xbp.  For 
xa(/j)  t°  actually  occur,  the  first  component  of  xa(ju)  must 
be  nonpositive.  That  is, 

xaM  =  {el)Tp{I-A)~lb<  0  (3) 

where  {el)T  =  (1  0  •••  0).  Similarly,  for  Xb(p)  to  actually 
occur,  one  needs 

x«,  (p)  =  {e1)Tp(I-B)~1b>  0.  (4) 


If  on  the  other  hand,  the  first  component  of  x,t  (p)  is  positive 
(the  first  component  of  xg(p)  is  nonpositive),  then  the  fixed 
point  is  called  a  virtual  fixed  point.  Virtual  fixed  points  are 
important  in  studying  the  dynamics  of  a  PWS  map  at  or 
near  the  border. 

Let  pa( X)  and  pb{X)  be  the  characteristic  polynomials 
of  A  and  B ,  respectively.  Then,  p,\(X)  =  detf/J  -  A)  and 
Pb(X )  =  det(Xl  —  B). 

The  fixed  points  can  be  written  as 


xa(p) 


(I  —  A)~1bp 
adj(7  —  A)bp 


det(7  — A) 

bA 

Pa(  if’ 


(5) 


and 


Xfi(p)  =  ( I  B )  'b/j 

adj  (I  —  B)b/j 
det  (I-B) 
bB 

Pb(  if’ 


(6) 


where  bA  =adj(/  —  A)b  and  bn  =adj(7  B)b.  It  can  be 
shown  that  b/\]  =  bBl  =:  b\  [4].  To  see  this,  recall  that  A 
and  B  differ  only  in  their  first  columns  and  adj  {I  —  A)  = 
(cof(I  —  A))T .  Thus,  the  first  row  of  adj (/  —  A )  is  equal  to 
the  first  row  of  adj(7  —  B),  which  implies  that  (e1)1 ad j(I  — 
A)b=  ( el)Tcidj(I~B)b=:b\ .  Thus,  the  first  component  of 
xa{b)  is  X/i,  (p)  =  jYvjP  anc*  t^le  ^rst  comPonent  of  xb(/j)  is 
xBl  (p)  =  pYjb-  F°r  t^le  fixed  point  xa(/j)  to  occur  for  p  <  0, 

it  is  required  that  xa1  (p)  <  0,  i.e.,  ®  ]>pT)  > 

Similarly,  for  the  fixed  point  xB{jP)  to  occur  for  ,u  >  0,  it 
is  required  that  xBi  (jj)  >  0,  i.e.,  >  0  >  0. 

Therefore,  a  necessary  and  sufficient  condition  to  have  a 
fixed  point  for  all  p  is  pA(l)ps(l)  >  0  which  is  assumed  to 
be  in  force  in  the  remainder  of  the  discussion. 

Next,  a  change  of  variables  is  performed  on  (2)  to 
simplify  the  analysis. 

Case  1):  p  <  0.  The  fixed  point  of  F  is  xA(/r)  =  (/— A)~lb/j. 
Changing  the  state  variable  in  (2)  to  z  =  x  —  xa  (p)  yields 
after  simplification 


z{k+  1) 


I  Az(k) ,  if  zi(k)  <  -xAl(p)  m 
\  Bz(k)+c/j ,  if  zi{k)  >  -xAi(b) 


where  c=  (B —A)(I  —A)~l)b.  In  the  new  coordinates,  z  =  0 
is  a  fixed  point  for  all  p  <  0.  (Note  that  the  border  Zborder  — 
{z :  zi  =  —  xAl  (p)}-  varies  as  a  function  of  p.)  Note  that  since 
B  and  A  differ  only  in  their  first  columns,  all  elements  of 
B  —  A  are  zero  except  for  the  first  column.  Thus,  c/u  =  (B  - 
A)(I—A)~lb^i=  (B-A)xa(Y)  =xa1(p)(B1  —A1),  where  the 
notation  A1  means  the  ;-th  column  of  the  matrix  A. 

Consider  the  quadratic  Lyapunov  function  candidate 

V  (z)  =  zTPz,  where  P  =  PT  >  0  (8) 


The  forward  difference  of  V  along  trajectories  of  (7)  is 
AV(z(k))  =  V(z(k  +  1)  —  V(z(k)).  There  are  two  cases: 
zi{k)  <  -xa,(p)  and  zi (k)  >  -xAl(ju). 

Case  1.1):  z i  ( k )  <  -xax  (p) 

A  VL{z{k))  =  V(z(*+l))-V(z(Jk)) 

=  (Az(k))T PAz(k)  -  z(k)T Pz(k) 

=  z(k)T(ATPA-P)z(k)  (9) 


Case  1.2):  zi(^)  >  -xax(ju) 

A  VR(z(k))  =  V(z(k+l))-V(z(k)) 

=  (Bz(k)  +  cju)TP(Bz(k)  +  cju)  - z(k)TPz(k) 

=  z(k)T(BTPB-P)z(k)  +  2jucTPBz(k)  +p2crPc 
=  z{k)T(BTPB  -  P)z{k)  +  2xAl  (p) (Bl  -A1)TPBz(k) 
+x2Al(Y(B'-Al)TP(B'-Al)  (10) 


Combining  (9)  and  (10)  yields 


A V(z(*)) 


A VL(z{k)),  if  zi(k)  <  -xAi(n)  nn 
AVr(z(^)),  if  zi  (k)  >  -xAl  0)  1 


From  (9)  and  (10),  a  necessary  condition  for  A V(z(k)) 
to  be  negative  definite  is  that  the  following  two  matrix 
inequalities  hold: 


AtPA-P  <  0,  (12) 

BtPB-P  <  0.  (13) 


Moreover,  the  following  claim,  which  asserts  sufficiency 
of  (12), (13)  for  negative  definiteness  of  A V{z{k))  is  stated 
and  proved. 


Claim:  (Sufficiency  of  LMIs  (12)-(13)  for  Decreasing  Lya¬ 
punov  Function) 

If  the  matrix  inequalities  (12)-(13)  are  satisfied  with  P  = 
PT  >  0,  then  A V(z(k))  given  by  (11)  is  negative  definite. 

Proof:  Assume  that  there  is  a  P  =  PT  >  0  such  that  (12)- 
(13)  are  satisfied.  Then,  AVz,  =  zT(ATPA  —  P)z  <  0  Vz  /  0. 
It  remains  to  show  that  A Vr  <  0.  Let  z  =  (zi,Z2)r,  where 
Zi  £  R  and  Z2  G  R"-1.  Note  that  AV  is  continuous  for 
all  z.  Continuity  of  AV  follows  from  the  continuity  of  V 
and  continuity  of  the  map  (7).  Since  A Vl  <  0  (AV/  =  0 
if  and  only  if  z  =  0)  and  AV  is  continuous  for  all  z,  it 
follows  that  A  Vr  <  0  at  the  border  {zi  =  — xAl(p)}  (since 

lim(zuZ2)^M(BU2)AVL  =  limLl,Z2)-(-^M,Z2)  AV«)-  11  re- 
mains  to  show  that  AVr  <  0  for  all  z  in  the  region  zi  > 

— xax  (p)  (note  that  — xAl  (p)  >  0).  Completing  the  squares 
in  (10)  allows  us  to  write  A Vr(z)  as  follows: 

AVr(z)  =  zT (BT PB  —  P)z  +  2xAx  {p)(Bl  —Al)TPBz 
+x2Al(Y(Bl-Al)TP(Bl-Al) 

=  (z  -  a )t(BtPB  -  P)  (z  -  a)  -  aT  ( BTPB  -P) a 
+xll(Y(Bl~Al)TP(B1-Al )  (14) 

where  a  =  -xAl(l-i)(BT PB  -  P)~lBT  P{Bl  -  A1).  Let  Nc 
R"  such  that  N  is  convex  and  contains  the  origin  (for 
example,  a  ball).  Since  the  fixed  point  xA(p)  is  close  to 
the  origin  for  small  p,  the  hyperplane  zi  =  —  xAl  (p)  slices 
the  neighborhood  N.  Consider  A Vr(z)  restricted  to  N.  The 
second  derivative  of  AV«(z)  with  respect  to  z  (i.e.,  its 
Hessian  matrix)  is  V2A Vr  =  2(BT PB  —  P)  <  0.  Thus,  AVr(z) 
is  strictly  concave  on  N,  i.e.,  for  every  z,  y  £  N,  and 
0  G  (0,1),  AVr(0z+  (1  -  9)y)  >  0A Vr(z)  +  (1  -  0)AV*(y). 
Note  that  AVr(0)  =xAi(p)(Z?1  —A,)TP(B1  —A1)  >  0.  Now, 
we  show  that  AVr  <  0  Vz  G  N  with  zi  >  — xAl  (p) .  By  way  of 
contradiction,  suppose  there  is  a  y  G  N,  with  y\  >  — xAi(p), 
such  that  AVtf(y)  >  0.  Since  A Vr(z)  is  strictly  concave, 
it  follows  that  AVr(z)  is  positive  along  the  line  segment 
connecting  0  and  v:  AVfl (0  •  0  +  (1  —  0)y)  >  0  Ays(0)  +(1  — 

>0 

0)AV/;(y)  >  0,  V0  G  (0,1).  But,  along  the  line  connecting 
>0 

Z  =  0  with  z  =  y,  there  is  a  point  z*  with  z*  =  — xAl  (p)  where 


AVR(z*)  <  0,  which  is  a  contradiction.  Thus,  A Vr(z)  <  0  for 
all  z  G  N  with  z  1  >  —  XAt  (p)  >  0.  ■ 

The  following  proposition  summarizes  the  results  so  far. 

Proposition  1:  The  forward  difference  of  V  =  zTPz,  with 
P  =  PT  >  0,  along  trajectories  of  (7)  with  p  <0  is  negative 
definite  (i.e.,  AT (z.)  <  0)  if  and  only  if  the  following  matrix 
inequalities  hold: 

Ar  PA  —  P  <  0,  (15) 

BtPB  P  <  0.  (16) 


Case  2):  p>0.  The  fixed  point  of  F  is  xB(p)  =  (/  —  B)~lbp. 
Changing  the  state  variable  in  (2)  to  z  =  x  —  xr{p)  yields 
after  simplification 


z(k+ 1) 


Az{k)  +  cp ,  if  zi  (k)  <  ~xBl  (ju)  nj) 
Bz(k),  if  zi(k)  >  -xBl(p) 


where  c  =  (A  —  B)  (/  —  B)  ~ 1  )b.  In  the  new  coordinates,  z  =  0 
is  a  fixed  point  for  all  p>0.  (Note  that  the  border  Zborder  = 
{z  ■  Zi  =  —  xBl(p)},  varies  as  a  function  of  p.)  Note  that 
since  B  and  A  differ  only  in  their  first  columns,  all  elements 
of  A  —  B  are  zero  except  for  the  first  column.  Thus,  cp  = 
(A -/?)(/ -B)~lbp  =  (A  —  B)xb(p)  =  xBl(p)(A1  -Bl). 

Consider  the  same  quadratic  Lyapunov  function  candi¬ 
date  as  in  (8)  above: 


V  (z)  =  zTPz ,  where  P  =  PT  >  0 


The  forward  difference  of  V  along  trajectories  of  (17) 
is  A V(z(k))  =  V(z(k+  1)  —  V(z(£)).  There  are  two  cases: 
Zi  ( k )  <  -xBl  ( p )  and  zi  ( k )  >  - xB{  (p)-  (Note  that  xB\ (p)  >  0 
from  (4).) 

Case  2.1):  zi  ( k )  <  —xBl  (p) 

A  VL(z(k))  =  v(z(k+\))-V(z{k)) 

=  (Az(k)  +  cp)T P(Az(k)  +  cp)  - z(k)T Pz(k) 

=  z(k)T  (Ar  PA  —  P)z(k)  +  2pcT  PAz(k)  +p2cTPc 
=  z(k)T(ATPA-P)z(k)  +  2xBl(p){A 1  -B1)TPAz(k) 
+x2Bl(p)(A1-Bl)TP(A1-Bl)  (18) 


Case  2.2):  zi  (k)  >  -xB]  ip) 

A  VR(z(k))  =  V(z(k+l))-V(z(k)) 

=  (Bz(k))TPBz{k)-z(k)TPz(k) 

=  z(k)T  (BtPB  —  P)z(k)  (19) 


Combining  (18)  and  (19)  yields 


A  V{z{k)) 


A VL{z{k)),  if  z(k)<-xBl{p) 

A VR(z(k)),  if  z(k)>-xBl(p)(M) 


Proposition  2:  (Necessary  and  Sufficient  Conditions  for 
Decreasing  Lyapunov  Function) 

The  forward  difference  of  V  =  zTPz,  with  P  =  PT  >0, 
along  trajectories  of  (17)  with  p  >  0  is  negative  definite  (i.e., 


AT  (z)  <  0)  if  and  only  if  the  following  matrix  inequalities 
hold: 

AtPA-P  <  0,  (21) 

BtPB  P  <  0.  (22) 


Proof:  Necessity  follows  from  (18)  and  (19),  and  the  proof 
for  sufficiency  is  similar  to  that  for  the  case  p  <  0  above.  ■ 

By  combining  Proposition  1  and  Proposition  2,  the  main 
result  of  this  paper  is  obtained. 

Proposition  3:  (Sufficient  Condition  for  Nonbifurcation 
with  Persistent  Stability  in  n-Dimensional  PWS  Maps ) 

The  PWS  map  (2)  has  a  globally  asymptotically  stable  fixed 
point  for  all  p  £  R.  if  there  is  a  P  =  PT  >  0  such  that 

AtPA-P  <  0, 

BtPB  P  <  0. 


Corollary  1:  If  at  p  =  0  the  origin  of  the  map  (2)  is 
quadratically  stable,  i.e.,  using  a  quadratic  Lyapunov  func¬ 
tion  V  =  xTPx,  with  P>  0,  then  the  fixed  point  depending  on 
p  on  both  sides  of  the  border  is  attracting  and  no  bifurcation 
occurs  from  the  origin  as  p  is  varied  through  zero. 

Below,  a  numerical  example  is  given  to  demonstrate  how 
the  Lyapunov-based  techniques  considered  in  the  previous 
section  can  be  used  in  the  stability  and  bifurcation  analysis. 


Example  1:  Consider  the  three-dimensinal  PWS  map 


x(k  + 1 ) 

f  Ax(k)  +  bp,  xi(k)<0 
}  Bx(k)  +  bp,  xi(A:)>0 

where 

(  0.4192 

0.3514 

0.3473  N 

A  =  1 

0.2840 

-0.2733 

-0.3107 

^  0.1852 

-0.2224 

-0.3974  ) 

(  -0.60 

0.3514 

0.3473  \ 

B  = 

0.56 

-0.2733 

-0.3107 

and  b  = 

^  -0.90 

-0.2224 

-0.3974  j 

The 

eigenvalues  of 

A  and  B  are  < 

(23) 


{0.5653,-0.7413,-0.0755} 
{0.0395,-0.6551  ±  j  0.4246}, 
both  A  and  B  are  Schur  stable 


o(A)  = 
and  = 

respectively.  Although 
matrices,  it  cannot  be 


concluded  that  no  bifurcation  for  (23)  occurs  at  p  =  0. 

A  common  quadratic  Lyapunov  function  V  =  xTPx,  with 
P  =  PT  >  0  that  satisfies  the  conditions  of  Proposition  3 
exists  for  this  example.  To  wit: 


/  1.6304  0.1559  -0.1313  \ 

P=  0.1559  1.3200  0.4436 

\  -0.1313  0.4436  1.3266  / 


is  obtained  using  the  MATLAB  LMI  toolbox.  Thus,  the 
PWS  map  (23)  has  a  unique  attracting  fixed  point  for  all  p 
(see  Figure  1). 


-0.2  — 
-0.2 


-0.1 


0  0.1  0.2 


Fig.  1.  Bifurcation  diagram  for  Example  1.  Each  solid  line  represents  a 
path  of  stable  fixed  points. 


IV.  Lyapunov-Based  Feedback  Control  Design 

In  this  section,  the  results  of  Section  III  are  used  in  the 
design  stabilizing  feedback  control  laws.  It  is  important 
to  emphasize  that  for  our  approach  to  apply,  the  control 
action  should  not  introduce  discontinuity  in  the  map.  This 
is  because,  as  discussed  in  the  introduction,  the  definition 
of  BCBs  requires  that  the  system  map  be  continuous  at 
the  border,  and  thus  our  results  on  nonbifurcation  with 
persistent  stability  also  apply  only  under  this  condition. 
Therefore,  to  maintain  continuity  of  the  map  after  control 
is  applied,  it  is  assumed  that  the  input  vectors  on  both  sides 
of  the  border  are  equal.  In  this  work,  the  input  vectors  are 
taken  to  be  equal  to  b  (the  derivative  of  the  map  with  respect 
to  the  bifurcation  parameter.) 

Simultaneous  feedback  control  is  considered  first,  fol¬ 
lowed  by  switched  feedback  control  design. 


A.  Simultaneous  Feedback  Control  Design 

In  this  control  method,  the  same  control  is  applied  on 
both  sides  of  the  border.  The  purpose  of  pursuing  stabilizing 
feedback  acting  on  both  sides  of  the  border  is  to  ensure 
robustness  with  respect  to  modeling  uncertainty.  Moreover, 
transformation  to  the  normal  form  is  not  required  when 
simultaneous  control  is  used.  All  that  is  needed  is  a  good 
estimate  of  the  Jacobian  matrices  on  both  sides  of  the 
border. 

Consider  the  closed-loop  system  using  static  linear  state 
feedback 


x(k  +  1) 
u{k) 


(  Ax(k)+bp  +  bu(k),  if  xi(k)<0 
\  Bx(k)  +  bp  +  bii(k) .  if  „n(k)>0 
gx(k)  (25) 


where  g  is  the  control  gain  (row)  vector. 

The  following  proposition  gives  stabilizability  condition 
for  the  border  collision  bifurcation  with  this  type  of  control 
policy. 


Proposition  4:  If  there  exists  a  P  =  PT  >  0,  and  a  feed¬ 
back  gain  (row)  vector  g  such  that 


P-(A  +  bg)TP(A  +  bg)  >  0  (26) 

P-(B  +  bg)TP(B  +  bg)  >  0  (27) 


then,  any  border  collision  bifurcation  that  occurs  in  the 
open-loop  system  (u  =  0)  of  (24)  can  be  eliminated  and 
persistent  stability  is  guaranteed  using  simultaneous  feed¬ 
back  (25).  Equivalently,  if  there  exists  a  Q  and  y  such  that 


(  Q  AQ  +  by 

\  ( AQ  +  by)T  Q 

(  Q  BQ  +  by 

\  ( BQ  +  by)T  Q 


>0, 

>0, 


(28) 

(29) 


then  any  border  collision  bifurcation  that  occurs  in  (24) 
can  be  eliminated  using  simultaneous  feedback  (25).  Here 
Q  =  P~l  and  the  feedback  gain  is  given  by  g  =  yP. 


Proof:  The  closed-loop  system  is  given  by 


x(k  +  1 )  = 


(A  +  bg)x(k)  +pb,  if  x\  (, k )  <  0 
(B  +  bg)x(k)  +pb,  if  x\  (k)  >  0 


(30) 


Using  Proposition  3,  a  sufficient  condition  to  eliminate  the 
BCB  is  the  existence  of  a  P  =  PT  >  0  such  that 


P  —  (A  +  bg)T  P  {A  +  bg)  >  0  (31) 

P-(B  +  bg)rP(B  +  bg)  >  0  (32) 


where  g  is  the  control  gain  to  be  chosen. 

Next,  inequalities  (31)-(32)  are  shown  to  be  equivalent 


to  (28)-(29)  using  the  Schur  complement  [25], 

[22], 

It 

straightforward  to  show  that 

P  —  (A  +  bg) T P  (A  +  bg) 

> 

0 

<=>  p-l-(A  +  bg)P~l(A  +  bg)T 

> 

o, 

and  P  —  (B  +  bg)  TP(B  +  bg) 

> 

0 

«=>  P-1-(B  +  bg)P~\B  +  bg)T 

> 

0. 

The  nonlinear  matrix  inequalities  above  are  transformed  into 
LMIs  using  the  Schur  complement  [25],  [22]: 


P-1  -(A  +  bg)p-\A  +  bg)T 

=  P-1  -{A  +  bg)p-xPP-\A  +  bg)T 

=  P'1  -  (AP-1  +bgP-')P(Ap-x  +bgP~l)T  >  0 

P*1  AP-x+by 
(. AP~l+by)T  P-1 


>0 


Similarly, 


P-{B  +  bg)TP(B  +  bg) 
P-1  BP  1  +  by  \ 

(BP-l+by)T  P-1  ) 


>  0 
>  0 


by  similar  reasoning.  ■ 

The  following  proposition  states  that  if  a  CQLF  exists  in 
one  coordinate  system,  another  CQLF  exists  in  a  different 
coordinate  system  arrived  at  using  a  simultaneous  similarity 
transformation  applied  to  both  A  and  B. 


Proposition  5:  (CQLF  and  Similarity  Transformations) 
Suppose  V  =  xTPx  (with  P  =  PT  >  0)  is  a  common  quadratic 
Lyapunov  function  for  both  of  the  matrices  A  and  B  (i.e., 
A'  PA  - P  <  0  and  BTPB -P<  0).  Then  V  =  xTPx  with  P  = 
(t^)tpt-1  =PT  >  0  is  a  common  quadratic  Lyapunov 
function  for  A  =  TAT~ 1  and  B  =  T BT  1  (i.e.  AT PA  —  P<  0 
and  Bt PB  —  P  <  0).  In  other  words,  if  a  CQLF  exists  in  one 
coordinate  system,  another  CQLF  exists  if  a  simultaneous 
change  of  coordinates  is  applied  to  both  A  and  B. 

Proof:  See  [16]. 

Remark  1:  The  switching  control  design  presented  above 
does  not  depend  on  the  border  separating  the  two  regions  of 
smooth  behavior.  Thus,  transformation  to  the  normal  form 
is  not  required  before  the  control  design. 


B.  Switched  Feedback  Control  Design 

Consider  the  closed-loop  system  using  static  piecewise 
linear  state  feedback 


f  Ax(k)  +  bp  +  bu(k) ,  if  X\ (k)  <  0 
{  Bx(k)  +  bp  +  bu(k),  if  xi(k)>0 


where 


gix(k),  x\ (k)  <  0 
g2x(k),  x\ (k)  >  0 


(33) 


(34) 


with  the  restriction  that  gi  and  g2  may  only  differ  in  their 
first  component,  i.e.,  g \ ,  =  g2j,  i  =  2,3,-  •  •  ,n.  This  condition 
is  imposed  to  maintain  continuity  along  the  border  {x :  x\  = 
0}. 


Proposition  6:  If  there  exists  a  P  =  PT  >0,  and  feedback 
gains  g i  and  g2  with  gy  =  g2i ,  i  =  2, 3,  •  •  • , n  such  that 


P-(A  +  bgi)TP(A  +  bgi)  >  0,  (35) 

P-(B  +  bg2)TP(B  +  bg2)  >  0,  (36) 


then  any  border  collision  bifurcation  that  occurs  in  the 
open-loop  system  ( u  =  0)  of  (33)  can  be  eliminated  using 
switching  feedback  (34).  Equivalently,  if  there  exists  a  Q, 
y i  and  ael  such  that 


(  Q 

\  m+byi)T 


(  Q  BQ  +  byA 
\m+byi)T  Q  ) 


-a 


AQ+Qby'  )>0, 

'  o  b(el)TQ\ 
yQelbT  0  ) 


(37) 
>  0,  (38) 


then  any  border  collision  bifurcation  that  occurs  in  (33)  can 
be  eliminated  using  switching  feedback  (34).  Here,  Q  =  P~l 
and  the  feedback  gains  are  given  by  gi  =  y\P  and  go  = 

g  t- a^r- 

Proof:  The  closed-loop  system  is  given  by 


f  (A  +  bgi)x(k)+pb,  if  x\ (k)  <  0 
\  (B  +  bg2)x(k)+pb,  if  xi(k)>0 


(39) 


Using  Proposition  3,  a  sufficient  condition  to  eliminate  the 
BCB  is  the  existence  of  a  P  =  PT  >  0  such  that 


P~(A  +  bgl)TP(A  +  bgl)  >  0  (40) 

P-(B  +  bg2)TP{B  +  bg2)  >  0  (41) 


where  g\,  g2  are  the  control  gains  to  be  chosen.  Inequali¬ 
ties  (40),(41)  are  equivalent  to 

Q  A  Q  +  by  i 

(AQ  +  byi)T  Q 

Q  BQ  +  by2 

C BQ  +  by2)T  Q 

respectively,  where  Q  =  P  ,  gi  =y\P  and  g2  =  y2P-  This 
equivalence  can  be  shown  using  similar  reasoning  as  that 
used  in  the  proof  of  Proposition  4. 

The  restriction  gy  =  g 2/,  i  =  2, 3,  •  •  •  ,n,  can  be  written  as 

g2  =  gi-a(el)T  (44) 


>  0  (42) 

>  0  (43) 


where  a  £  R.  Therefore, 

yi-yi  =  gxQ-giQ 
=  (gi-gi)Q 

=  a  (e1)TQ  (45) 

Substituting  y2  =  }’i  —  a (e1)TQ  in  (43)  yields  (38).  This 
completes  the  proof.  ■ 

Note  that  if  a  =  0  in  (38),  then  the  switching  feedback 
control  (34)  becomes  simultaneous  control. 

Remark  2:  We  remark  that  switching  control  design 
(with  no  restriction  on  feedback  gains)  was  used  in  [22] 
for  stabilization  of  the  origin  of  discrete  time  switching 
systems.  No  bifurcation  control  was  considered  in  the 
referenced  work. 


C.  Numerical  Examples 

Below,  numerical  examples  that  demonstrate  the  pro¬ 
posed  feedback  control  methods  are  given. 


Example  2:  (Fixed  point  attractor  bifurcating  to  instan¬ 
taneous  chaos) 

Consider  the  three  dimensional  PWS  map 


r(k+l)  = 


Ax(k)  +  bp,  x\  (, k )  <  0 
Bx(k)  +  bp,  x\  ( k )  >  0 


(46) 


where 


A  = 


B  = 


0.0334 

—0.4588 

0.0474 

0.8384 

-0.8180 

0.6602 


1.7874 

-0.4430 

-0.0416 

1.7874 

-0.4430 

-0.0416 

of 


-0.1705  \ 

-0.8282  , 

0.8000  ) 

-0.1705  \ 

-0.8282  and  b  = 
0.8000  ) 


The  eigenvalues  of  A  and  B  are  a  (A)  = 
{0.766,-0.1878  ±  j  0.8389}  and  a(B)  = 
{-0.1157,0.6555  ±  j  1.0987},  respectively.  Note  that 
A  is  Schur  stable,  but  B  is  unstable.  Simulation  results 
show  that  (46)  undergoes  a  border  collision  bifurcation 
from  a  fixed  point  attractor  to  instantaneous  chaos  at  p  =  0 
(see  Figure  2). 


0.2 


0.2 


-0.2 - 1 - 1 - ' - 

-0.2  -0.1  0  0.1  0.2 

Fig.  2.  Bifurcation  diagram  for  Example  2.  The  solid  line  represents  a 
path  of  stable  fixed  points  and  the  shaded  region  represents  a  one  piece 
chaotic  attractor  growing  out  of  the  fixed  point  at  /j  —  0. 


-0.2 - 1 - 1 - 1 - 1 
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Fig.  3.  Bifurcation  diagram  for  Example  2  with  simultaneous  feedback 
control  u(k)  —  gx(k).  The  solid  lines  represent  a  path  of  stable  fixed  points. 


Feedback  control  design:  Using  the  results  of  Proposition  4, 
a  symmetric  and  positive  definite  matrix  Q  and  a  feedback 
control  gain  vector  g  that  satisfy  the  LMIs  (28)-(29)  are 
sought.  A  solution  to  (28)-(29)  is  obtained  using  the  MAT- 
LAB  LMI  toolbox.  To  wit: 


/  0.4753 

-0.0428 

-0.1694 

Q 

=  -0.0428 

0.8821 

-0.1647 

\  -0.1694 

-0.1647 

0.5041 

y 

=  (  -0.1601  - 

-1.4937 

0.3356  ), 

8 

II 

V! 

to 

1 

=  (  -0.5193  - 

-1.7324 

-0.0747  ) 

closed-loop  matrices  are 

given  by 

Ac 

=  A +bg 

/  -0.4859 

0.0550 

-0.2452 

=  -0.4588 

-0.4430 

-0.8282 

\  0.0474 

-0.0416 

0.8000 

Be 

=  B  +  bg 

/  0.3191 

0.0550 

-0.2452 

=  -0.8180 

-0.4430 

-0.8282 

\  0.6602 

-0.0416 

0.8000 

Their  eigenvalues  are:  o(Ac.)  =  {0.8141,-0.4715  ± 
j  0.1409}  and  a(Bc)  =  {-0.4507, 0.5634 ±j  0.3498}.  The 
bifurcation  diagram  of  the  closed  loop  system  is  depicted 
in  Figure  3. 

Example  3:  (Saddle-node  border  collision  bifurcation ) 
Consider  the  three  dimensional  PWS  map 


/  3.3000 

-0.2280 

-0.9385  \  , 

f  1  \ 

B  =  1 

-0.6299 

-0.0029 

0.9191  j  and  b  =  j 

° 

l  0.3705 

-0.5107 

0.5553  /  ' 

)  0  / 

The  eigenvalues  of  A  and  B  are  a(A)  =  {—0.2921,0.4397  ± 
j  0.3470}  and  a(B)  =  {3.1739,0.3392  ±j  0.4756},  re¬ 
spectively.  Note  that  A  is  Schur  stable,  but  B  is  unstable. 
Simulation  results  show  that  (47)  undergoes  a  saddle  node 
border  collision  bifurcation  where  a  stable  and  an  unstable 
fixed  point  collide  and  disappear  as  p  is  increased  through 
zero  (see  Figure  4). 


Fig.  4.  Bifurcation  diagram  for  Example  3  without  control.  The  solid  line 
represents  a  path  of  stable  fixed  points  whereas  the  dashed  line  represents 
a  path  of  unstable  fixed  points. 


x(k  + 1 ) 


f  Ax(k)+bp1  ri(k)<0 
{  Bx(k)  +  bp,  ,*'i(k)>0 


where 

/ 

(  0.0350 

-0.2280 

-0.9385 

A  = 

-0.3123 

-0.0029 

0.9191 

1 

l  -0.3825 

-0.5107 

0.5553 

(471 

v  ’  Feedback  control  design:  Simultaneous  stabilizing  feedback 
control  based  on  Proposition  4  does  not  exist  for  this 
example.  Therefore,  a  stabilizing  switched  feedback  control 
using  Proposition  6  is  sought.  Using  the  LMI  toolbox  in 
MATLAB,  a  symmetric  and  positive  definite  matrix  Q  and 
a  feedback  control  gain  vectors  gi  and  g2  that  satisfy  the 


LMIs  (37)-(38)  are  obtained: 


a 

=  3.0972 

/  25.3606 

4.5507 

7.9810 

Q 

=  4.5507 

43.0961 

9.8713 

\  7.9810 

9.8713 

30.8840 

y  i 

=  (  5.7709 

14.8260 

34.4887  ) 

7 

Oj 

£ 

II 

=  (  -0.1436 

0.1024 

1.1211  ), 

82 

=  g\-u{el)T 

=  (  -3.2408 

0.1024 

1.1211  ). 

The  closed-loop  matrices  are  given  by 

Ac  =  A  +  bgi 


f  -0.1086 

-0.1256 

0.1826 

=  -0.3123 

-0.0029 

0.9191 

\  -0.3825 

-0.5107 

0.5553 

=  B  +  bg2 

/  0.0592 

-0.1256 

0.1826 

=  -0.6299 

-0.0029 

0.9191 

\  0.3705 

-0.5107 

0.5553 

Their  eigenvalues  are:  o(Af)  =  {0.0011,0.2213±  j  0.6236} 
and  o(Bc)  =  {-0.0002, 0.3059  ±j  0.5102}.  The  bifurcation 
diagram  of  the  closed-loop  system  is  similar  to  Figure  3. 

V.  Concluding  Remarks 

Lyapunov-based  analysis  of  piecewise  smooth  discrete¬ 
time  systems  that  undergo  border  collision  bifurcations 
has  been  considered.  One  of  the  main  contributions  of 
this  paper  is  that  a  sufficient  condition  for  nonbifurcation 
with  persistent  stability  in  PWS  maps  of  dimension  n  that 
depend  on  a  parameter  was  derived.  This  condition  has 
been  used  in  the  design  of  stabilizing  feedback  control 
laws  to  eliminate  border  collision  bifurcations  in  PWS  maps 
and  produce  desirable  behavior.  Numerical  examples  were 
given  to  demonstrate  the  efficacy  of  the  proposed  control 
techniques. 
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